Millisecond pulsars with r-modes 
as steady gravitational radiators 
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Millisecond pulsars (MSPs) are generally agreed to originate in low-mass X-ray binaries (LMXBs), 
in which mass transfer onto the neutron stars spins them up to their observed, fast rotation. The 
lack of MSPs and LMXBs rotating near break-up and the similar rotation periods of several LMXBs 
have been attributed to the accretion torque being balanced, at fast rotation, by gravitational radi- 
ation, perhaps caused by an oscillation mode made unstable through the so-called Chandrasekhar- 
Friedman-Schutz mechanism. Recently, Wagoner has argued that internal dissipation through weak 
interaction processes involving A and E~ hyperons may cause LMXBs to evolve into a quasi-steady 
state, in which the neutron star has a nearly constant rotation rate, temperature, and mode ampli- 
tude. We take this hypothesis one step further, showing that MSPs descending from these LMXBs 
spend a long time in a similar state, in which a low-amplitude r-mode turns them into extremely 
steady sources of both gravitational waves and thermal X-rays, while they spin down due to a combi- 
nation of gravitational radiation and the standard magnetic torque. Observed MSP braking torques 
already place meaningful constraints on the allowed gravitational wave amplitudes and dissipation 
mechanisms. 

PACS numbers: 04.30.Db, 97.60.Gb, 97.60.Jd, 97.80.Jp 



The Chandrasekhar-Friedman-Schutz instability Q, 0, 
3] can make retrograde r-modes (inertial modes due to 
the Coriolis force) on rapidly rotating neutron stars grow 
and emit gravitational waves at the expense of the stel- 
lar rotational energy and angular momentum 0, • R 
is choked |(| at very high temperatures by bulk viscos- 
ity due to non-equilibrium weak interactions such as the 
Urea processes (n — * p+e+P e andp+e — * n+i> e , where n, 
p, e, v e , and v e denote neutrons, protons, electrons, elec- 
tron antincutrinos, and electron neutrinos, respectively, 
or analogous processes involving hyperons), and at very 
low temperatures by standard shear viscosity. Additional 
damping could be provided at intermediate temperatures 
by hyperon bulk viscosity 0, @ , caused by the processes 
A + n n + n and E~ + p <-> n + n. If strong enough, 
this can split the instability region on the temperature- 
frequency plane into two separate windows, as illustrated 
in Fig. 1. 

In the standard "recycling" model, millisecond pulsars 
(MSPs) reach their high rotation rates due to the accre- 
tion torque in a low-mass X-ray binary system (LMXB) 
0. For this to happen, some mechanism must suppress 
the CFS instability enough to prevent the gravitational 
wave torque from being larger than the accretion torque. 
The simplest scenario [ill requires a shear viscosity sub- 
stantially higher than that used in Fig. 1, so that the 
instability boundary at the equilibrium temperature T acc 
reached during the accretion process 11] lies higher than 
the fastest rotation frequencies observed in MSPs, ~ 650 
Hz. This is not unlikely, since a stronger effective shear 
viscosity might be obtained with a thicker crust [l^ , su- 
perfluids [13j, or a magnetic field |l4|. An instability 
boundary just above the frequency of the fastest pulsars 



might even explain why there are no MSPs approaching 
the break-up frequency, but it would tend to predict an 
accumulation of MSPs near this boundary, rather than 
the observed scarcity of MSPs with periods below 2 ms 

Wagoner has shown that fast rotations periods 
can be reached even if there is a substantial unstable 
region at T acc . In this scenario, the low-temperature in- 
stability window is reached from below as the neutron 
star is spun up by accretion (see Fig. 2(a)). A stable 
balance between cooling through neutrino emission and 
heating through nuclear reactions undergone by the ac- 
creted matter in the deep crust (which release ~ 1 MeV 
per baryon |ll|) keeps the stellar interior warm, at a 
temperature T acc ~ ]_o 7 - 5-8 - 8 K that depends on the al- 
lowed neutrino emission processes (e.g., direct vs. mod- 
ified Urea). As the star moves into the unstable region, 
its most unstable r-mode 0| (with spherical-harmonic 
indices I = m = 2) is excited, producing viscous dis- 
sipation. The resulting heat release moves the star to 
the high-temperature boundary of the instability window 
[3G3, set by 



0.051^ Hz s-i = 0.13/^ k 2 Hz T| s 
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The left-hand side gives the r-mode driving rate through 
gravitational radiation reaction 6] for a "fiducial" neu- 
tron star with a polytropic (n = 1) equation of state, 
mass M = 1.4 Mq, and radius R — 12.5 km, in terms 
of the rotation frequency in kHz, i-'kHz- The right-hand 
side is the damping rate, dominated by hyperonic pro- 
cesses, whose large uncertainty (including the reduction 
due to superfluid effects) is parameterized by the dimen- 
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FIG. 1: Instability boundaries for the I = m = 2 r-mode. 
The curves all use the same shear viscosity, determined by 
Levin & Ushomirsky |2 to be active at the boundary be- 
tween the liquid core and a thin, elastic crust, but different 
assumptions about the bulk viscosity in the core. The thin 
dotted curve considers only modified Urea processes, whereas 
the other three also include direct Urea processes and the hy- 
peron bulk viscosities proposed by Jones Q (thin solid line) 
and by Lindblom & Owen Q (thick solid line), in both cases 
ignoring superfluid effects, and by Lindblom & Owen under 
the effects of hyperon superfluidity with a uniform, high crit- 
ical temperature, T c ~ 5 x 10 9 K (thick dotted line). In each 
case, the unstable region lies above the curve. 
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FIG. 2: Evolutionary tracks of LMXBs and MSPs. Panel (a) 
shows the evolution of a neutron star in a LMXB being spun 
up by accretion|lq|. and panel (b) that of a MSP being spun 
down by gravitational radiation and/or a magnetic torque 
(solid lines) . The instability boundary (dashed) considers the 
same shear viscosity as in Fig. 1, and the hyperon bulk vis- 
cosity normalized to fit the parameters of PSR B1957+20 (see 
text). 



sionless parameter fa. This parameter is unity for the 
preferred bulk viscosity of Lindblom and Owen Q with 
no superfluid effects, acting in a hyperonic core of radius 
r = R/2, whereas fh = 0.16 corresponds to the bulk vis- 
cosity previously proposed by Jones Q. The resulting 
interior temperature, in units of 10 8 K, turns out to be 
T 8 = 0.63/- 1/2 ^ Hz . 



Once the temperature is high enough for accretion 
heating to be negligible, the dimensionless mode ampli- 
tude |6j a is set by the balance of the hyperon bulk viscos- 
ity heating and the neutrino emission (the latter assumed 
to be dominated by direct Urea processes |l8j|). 

3.6 x 10 50 f h Tg4 Hz a 2 erg s" 1 = 1.2 x 10 39 /t/T 8 6 erg s~\ 

(2) 

where we again parameterize our ignorance by a dimen- 
sionless factor, fu, which is unity when the threshold 
for nucleonic direct Urea processes (n «-> p + e~) is 
reached at half the stellar radius. It can take different 
(possibly temperature-dependent) values if the thresh- 
old radius is different, if only hyperonic direct Urea pro- 
cesses are allowed |lj|, or if superfluidity reduces the re- 
action rates [2(| • The resulting amplitude is quite small, 
a = 7.5 x 10~ 7 {fu / fh) so the linear approxima- 

tion used in the expressions above are amply justified, 
and it appears unlikely that the mode could be satu- 
rated at lower amplitudes by any other processes such as 
turbulent viscosity [2l| or mode coupling [22j |. 

With the temperature and mode amplitude at any 
given rotation rate set by eqs. (1) and (2), the accre- 
tion torque, My/GMR = 1.5 x 10 34 (M/M E dd) dyne cm 
(where M is the mass accretion rate and MEdd = 
1.9 x 10~ 8 M Q yr _1 is the standard "Eddington limit" 
gradually drives the star "up" towards faster ro- 
tation. It may eventually settle in a quasi-steady state, 
in which the accretion torque is balanced by the angu- 
lar momentum loss due to gravitational radiation [16j, 
6.8 x 10 46 ^ Hz a 2 dyne cm, at an equilibrium period 

/ ' \ 1/11 

Peq= ^ =Ll U^"J ms ' ( ) 

which is only very weakly dependent on the uncertain 
parameters. In this state, the neutron star is a source 
of low-amplitude gravitational waves, whose steadiness 
depends on that of the accretion rate. 

It is worth noting that eq. (3) and the expression above 
for the accretion torque are only valid if the star's rota- 
tion is much slower than the mass-shedding limit and 
the accretion disk extends down to the stellar surface, 
not being constrained by the star's magnetic field. We 
also point out that the possibility of accretion torques in 
LMXBs being balanced by gravitational wave emission, 
in somewhat different scenarios, had previously been sug- 
gested Ulllij. 

In what follows, we analyze what happens in Wagoner's 
scenario once the accretion stops and the star turns on 
as a (millisecond) radio pulsar. It slowly slides down 
the instability boundary as it loses angular momentum 
through the combined effects of gravitational radiation 
and standard magnetic braking (Fig. 2(b)). Its heat- 
ing and cooling time scales, as well as the growth and 
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damping times for the r-modes, are much shorter than 
the spin-down time, and therefore the star will remain 
in the quasi-steady state described by equations (1) and 
(2). Since no external perturbations are present, the spin- 
down is extremely smooth and predictable, as observed 
in MSPs j2|| . For a given neutron star model (which 
sets the "fudge factors" fu and fh), the rotation pe- 
riod, P (observed as the periodicity of the radio pulses), 
uniquely determines the gravitational-wave contribution 
to the spin-down rate 



iW = 5.3X10-$ 

Jh 



lms 
~P~ 



(4) 



We note that, if the particles in the neutron star core are 
superfluid, the factors fu and fh will depend on temper- 
ature, and therefore implicitly on P, which would mod- 
ify the dependence of Pgw from the proportionality to 
P~ 9 implied by eq. (4). However, for typical cases that 
we evaluated numerically, a strongly decreasing function 
Pgw(P) is still obtained. 

Of course, for periods P of actually observed MSPs 
in the range of applicability of this model, Pgw (P) can 
be at most as high as the observed P. The most con- 
straining case is the so-called "black widow" ,27], PSR 
B1957+20, because of its short period. P = 1.60 ms, and 
its small intrinsic period derivative [2j|, P — 1.2 x 10 -20 . 
For this pulsar, we require fj*/fu > 6.4, i.e., a some- 
what larger bulk viscosity and/or a substantially lower 
neutrino emissivity than in our favored model. Neither 
of these would be too surprising, in view of their strong 
dependence on the uncertain state of very dense matter. 
The presence of superfluid energy gaps for the particles in 
the neutron star core could substantially reduce the phase 
space and therefore the rates of all reactions, thereby de- 
creasing both the Urea cooling rate [2(J and the hyperon 
bulk viscosity H, Eil, potentially down to fu,fh <C 1. 
Since the reactions determining fh involve more poten- 
tially superfluid particles (both hyperons and nucleons) 
than the Urea processes (just nucleons), it is likely that 
fh < fu, making it difficult to satisfy the "black widow" 
constraint. However, fu may be reduced if, instead of nu- 
cleonic direct Urea processes, only hyperonic direct Urea 
reactions are allowed, which would set fu (without su- 
perfluid reduction) at ~ 0.2 0. The latter processes are 
almost certainly allowed as long as hyperons exist in the 
star, but again may be reduced (perhaps more strongly 
than the neutrino- less processes contributing to fh) if 
the hyperons become superfluid at high temperatures. 
Therefore, instead of ruling out the model, the observa- 
tions may be used to constrain the internal properties of 
neutron stars once the model is validated in other ways. 

Of course, the values of fh and fu are most likely not 
identical for all MSPs, as these probably have different 
masses, and therefore different mean densities and differ- 
ent fractions of their interior in which the relevant pro- 



cesses take place. Assuming, for simplicity, that the con- 
dition ffjfu > 6.4 does apply to all MSPs and LMXBs, 
we can in principle constrain the LMXB equilibrium pe- 
riod to P cq < 0.9(M E dd/M) 1/n ms. Strictly speaking, 
this period is beyond the mass-shedding limit for our 
adopted neutron star model, and should therefore not be 
taken literally. It shows that, in order to be consistent 
with MSP observations, the r-mode instability cannot 
limit the rotation of LMXBs to much below the maxi- 
mum allowed equilibrium rotation for any neutron star, 
unless magnetic stresses prevent the accretion disk from 
reaching the surface, reducing the accretion torque and 
increasing the equilibrium period. A similar conclusion 
has been reached before |3Jj], based on the low luminos- 
ity of LMXBs in quiescence, compared to the expectation 
from internal dissipation. 

The spin-down of a newly born MSP, initially spinning 
at a period Pq near the mass-shedding limit, is likely 
to be dominated by the gravitational radiation torque. 
The time it takes to spin down to a slower period P 
is t « 3.9 x 10 6 (/3/6.4/ l/ )[(P/ms) 10 - (P /ms) 10 ] yr, 
substantially shorter (and much more period-dependent) 
than for magnetic braking at the inferred magnetic field 
strengths, ~ 10 8 ~ 9 G. This makes it much less likely to 
find such extremely rapidly spinning MSPs, perhaps in 
this way explaining the observed absence of these ob- 
jects 15] . An independent test of this scenario would 
be a measurement of the braking index, n = vv/v 2 |23l |. 
which reaches a huge value (n = 11, if fh/fu is con- 
stant) in the regime where the torque is dominated by 
gravitational- wave emission (compared to n = 3 for pure 
magnetic dipole braking). Unfortunately, measurements 
of braking indices in MSPs are out of reach. 

In this state, the effective surface temperature of the 
MSP can be inferred from its interior temperature [3l| 
to be T s w 8 x 10 5 / h ; - 28 (P/ms)- 11 K, not far below 
the upper limits obtained from ROSAT observations of 
MSPs ,32} , and perhaps within the reach of more careful 
determinations based on XMM-Newton spectra. 

Perhaps most interestingly, the MSP will radiate grav- 
itational waves at the r-mode frequency. In the limit 
of slow rotation and weak gravity, this frequency is 4/3 
times the rotation frequency |33J, or w 850 Hz for the 
fastest pulsars. A more precise determination, including 
centrifugal and relativistic corrections, should be possi- 
ble. The angle-averaged amplitude at a distance D is 
M 



7 x 10 



-28 



6.4/i 
fh 



1/2 



1.6ms\ 1.5kpc 



P 



D 



(5) 



where the reference numbers correspond to the "black 
widow" pulsar, the most favorable known so far. Ac- 
cording to recent sensitivity curves |35| . this signal is 
almost within reach of the advanced LIGO with signal 
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recycling, tuned at the appropriate frequency and inte- 
grating for 1/3 yr. Therefore, future gravitational-wave 
observatories may well give us information about weak 
interaction processes in superdense matter. 

A recent study j3|| has found that "strange stars", 
composed of "deconfined" u, d, and s quarks, may have 
a time evolution qualitatively similar to that later found 
by Wagoner 0] for neutron stars in LMXBs and here 
for neutron stars as MSPs. In the corresponding strange 
stars, the bulk viscosity is dominated by the process 
u + d <-> s +u, the quark analog for the hyperonic pro- 
cesses considered in the present discussion. In that study, 
the reheating due to bulk viscous dissipation is ignored, 
and therefore the evolution time scale is set by passive 
cooling, which is much faster and thus more violent than 
in the more realistic case discussed here. In view of the 
present theoretical uncertainties and lack of observational 
constraints, LMXBs and MSPs may well be strange stars 
rather than neutron stars. Although the numerical de- 
tails are different (and even more uncertain), "strange 
MSPs" may follow a similar evolutionary path and be 
subject to observational constraints analogous to those 
discussed here for neutron stars. 

We conclude that, if there is a substantial instability 
window at low temperatures, neutron stars (or strange 
stars) in LMXBs can generally spin up to MSP rota- 
tion rates by the mechanism suggested by Wagoner (lij . 
The gravitational wave torque will still be active in their 
MSP phase, until they reach the bottom of the instability 
window, which some of the observed MSPs may not yet 
have done. Thus, these MSPs will be sources of gravita- 
tional waves and thermal X-rays. Even if born rotating 
near break-up, MSPs will spin down to near the observed 
periods much more quickly than in the magnetic braking 
model, perhaps explaining the scarcity of very fast MSPs. 
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